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ABSTRACT 

We study formation and evolution of bar-disk systems in fully self-consistent cosmological simula- 
tions of galaxy formation in the ACDM WMAP3 Universe. In a representative model we find that 
the first generation of bars form in response to the asymmetric dark matter (DM) distribution (i.e., 
DM filament) and quickly decay. Subsequent bar generations form and are destroyed during the ma- 
jor merger epoch permeated by interactions with a DM substructure (subhalos). A long-lived bar is 
triggered by a tide from a subhalo and survives for ~ 10 Gyr. The evolution of this bar is followed 
during the subsequent numerous minor mergers and interactions with the substructure. Together with 
intrinsic factors, these interactions largely determine the stellar bar evolution. The bar strength and 
its pattern speed anticorrelate, except during interactions and when the secondary (nuclear) bar is 
present. For about 5 Gyr bar pattern speed increases substantially despite the loss of angular mo- 
mentum to stars and cuspy DM halo. We analyze the evolution of stellar populations in the bar-disk 
and relate them to the underlying dynamics. While the bar is made mainly of an intermediate age, 
~ 5 — 6 Gyr, disk stars at 2; = 0, a secondary nuclear bar which surfaces at z ~ 0.1 is made of younger, 
~ 1 — 3 Gyr stars. 

Subject headings: cosmology: dark matter — galaxies: evolution — galaxies: formation — galaxies: 
halos — galaxies: interactions — galaxies: kinematics and dynamics 



1. INTRODUCTION AND NUMERICS 

Within the framework of structure formation in the 
universe, a hierarchy of dark matter (DM) masses form, 
while baryons assemble in their midst as galactic disks 
(e.g.. White & Rees 1978). Disk evolution is expected 
to be influenced heavily by mergers and interactions in 
various ways. In this Letter, we focus on the specific 
issue of a tidal triggering of stellar bars embedded in 
a disk and surrounded by a DM halo with a substruc- 
ture, within the context of cosmological evolution in the 
WMAP3 universe (see also Dubinski et al. 2008). We an- 
alyze the dynamical aspects of this evolution and follow 
the prevailing stellar populations in the bar-disk system. 

Recent efforts to understand galaxy formation have 
been spearheaded by the work on pure DM halo for- 
mation (e.g., Diemand et al. 2007) and associated disk 
growth (e.g., Sommer-Larson et al. 2003; Governato et 
al. 2004, 2007; Heller et al. 2007a,b). Addition of bary- 
onic component (s) to this modeling has met with difficul- 
ties, partly due to a numerical resolution, star formation 
(SF), energy feedback, and other processes. Dissipation 
is relatively slow in disks, but can be accelerated by any 
substantial departure from axial symmetry due to grav- 
itational torques — a nonlocal viscosity, surpassing the 
conventional viscosity by orders of magnitude. Torques 
facilitate the redistribution of mass, angular momentum 
and energy between the disk baryons and the slowly tum- 
bling DM halo. Two main sources can serve as triggers 
of gravitational torques onto the disk — an asymmetric 

^ Department of Physics and Astronomy, University of Ken- 
tucky, Lexington, KY 40506-0055, USA 

2 JILA, University of Colorado, Boulder, CO 80309, USA 

^ Department of Physics, Georgia Southern University, States- 
boro, GA 30460, USA 

^ Racah Institute of Physics, Hebrew University; Jerusalem 
91904, Israel 



mass distribution in the DM halo, and its non- uniformity, 
i.e., substructure, in DM and baryons. Both tidal inter- 
actions (e.g., Byrd et al. 1986; Noguchi 1987; Gerin et al. 
1990) and DM halo asymmetry (Heller et al. 2007a,b) 
are likely to trigger stellar bars. About 1/3 of the large- 
scale bars host additional nuclear bars (Laine et al. 2002; 
Erwin & Sparke 2002). 

Observations and numerical simulations point to an 
intricate relation between the stellar populations and the 
underlying dynamics of bars (e.g., Martin & Friedli 1997; 
Knapen et al. 1995a, b; Jogee et al. 2002a, b; Zurita & 
Perez 2008). For bars in a cosmological context, this issue 
becomes even more relevant, because the disk is open to 
accretion of cold gas and energy input from interactions. 
Here we attempt to answer some questions about the 
long-term evolution of bar-disk populations and relate it 
to the dynamical history of the system. 

So far, nearly all modeling of bar evolution has been 
limited to isolated systems in a stationary state, with the 
exception of simulations of tidal interactions in controlled 
experiments (e.g., Gauthier et al. 2006). We test the bar 
origin and evolution in the cosmological setting without 
fine-tuning and minimazing the prior assumptions. 

Numerical simulations have been performed using the 
FTM-4.5 hybrid 7V-body/SPH code (e.g.. Heller & Shlos- 
man 1994; Heller et al. 2007b; Romano-Diaz et al. 
2008a) using physical coordinates. The total number 
of DM particles is 2.2 x 10^ and the SPH particles 
4 X 10^. The gravity is computed using the falcON rou- 
tine (Dehnen 2002) which scales as 0(N). The gravita- 
tional softening is 500 pc, for DM, stars and gas. We 
assume the ACDM cosmology with WMAP3 parame- 
ters, Qrn = 0.24, Qa = 0.76 and h = 0.73, where h is 
the Hubble constant in units of 100 km s~^Mpc~^. The 
variance = 0.76 of the density field convolved with the 
top hat window of radius Sh~^ Mpc is used to normalize 
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the power spectrum. The SF algorithm is described in 
Heher et al. (2007b). We note only that multiple gen- 
erations of stars can form from a single SPH particle. 
The energy and momentum feedback into the ISM are 
implemented, and the associated parameters have the 
following values (Heller et al.): energy thermalization 
esF = 0.3, cloud gravitational collapse ceff = 1, and self- 
gravity fudge- factor acrit = 0.5. 

Initial conditions generated here are those of Romano- 
Diaz et al. (2008a): we use the Constrained Realiza- 
tions (CR) method (Hoffman & Ribak 1991) within a re- 
stricted box size of 8h~^ Mpc and a sphere of 5h~^ Mpc 
is carved out and evolved. A large-scale filament runs 
across the computational sphere and is banana-shaped. 
The constructed Gaussian field is required to obey a 
set of constraints of arbitrary amplitudes and positions 
(more about this method in Romano-Diaz et al. 2006, 
2007). Two constraints were imposed on the density 
field, first — that the linear field Gaussian smoothed with 
a kernel of 1.0 x 10^^ /i^^M©, has an over-density of ^ = 3 
at the origin (2.5cr perturbation, where cr^ is the variance 
of the appropriately smoothed field). It was imposed on 
a 256^ grid and predicted to collapse at Zc ^ 1.33 based 
on the top-hat model. This perturbation is embedded 
in a region (2nd constraint) corresponding to a mass of 
5 X 10^^ h~^MQ in which the over- density is zero, i.e., 
the unperturbed universe. The random component in 
CRs favors formation of similar structures which leads 
to major mergers. The DM mass inside the computa- 
tional sphere is~6.1xl0-^^ h~^MQ. We have randomly 
replaced 1/6 of DM particles by equal mass SPH parti- 
cles. Therefore, Q^n is not affected. 

2. EVOLUTION AND STELLAR POPULATIONS 

The prime halo evolves through a series of major merg- 
ers which end by z ~ 1.5, while accretion of subhalos and 
smooth accretion continue to the present, mainly along 
the filament. The number of subhalos varies with time, 
and at its peak is about few hundred, for subhalos in 
excess of 10^ Mq. The disk is recognizable from ^ 8, 
associated with a strong SF. It grows from inside out and 
appears strongly non-axisymmetric and barred along the 
filament. The first merger of nearly equal masses is gas- 
rich, coplanar and prograde. The disk reforms immedi- 
ately, although the bar does not {z ^ 5). Next merger, 
at 2: ~ 4, is again prograde and excites a bar which is de- 
stroyed shortly. Subsequent minor mergers drive a strong 
spiral structure, as the disk grows to ~ 10 kpc, and push 
the gas inwards in a 'shepherding' mode. The SF is vig- 
orous in the high surface density gas. At z ~ 2.6, the 
tide from a subhalo (prograde and inclined penetrating 
minor merger) triggers a strong gas-rich bar (Fig. 1) fol- 
lowed by (minor) merger activity till z = 0. A direct hit 
by another disk (in a host subhalo), on ~ 40°-inclined 
prograde trajectory, weakens the bar abruptly at 2: ~ 2.2 
and shortens it. A retrograde encounter at z ~ 1.4 has 
no dramatic effect on the bar. As a result of these in- 
teractions the gas-rich disk shrinks gradually. Overall, 
the disk axis is oriented along the halo minor axis, but 
this orientation experiences strong departures during the 
major merger epoch. The DM halo forms an isother- 
mal cusp, R~^, within 15 kpc which survives beyond 
z 1 and is ultimately washed out, leaving a flat core of 
- 2 - 3 kpc (Romano-Diaz et al. 2008b). 



Bar dynamics is displayed in Fig. 2. The pattern speed, 
l^b, of the initial barlike response to the DM stagnates 
initially, then accelerates. The bar amplitude, A2 of the 
m = 2 Fourier mode, decays rapidly. The oval pertur- 
bation in the stellar disk remains but is weak and noisy, 
and therefore omitted here until the tidal triggering at 
z ~ 2.6. Subsequent evolution depends both on interac- 
tions with the substructure and on intrinsic factors. 

For the newly triggered bar, l^b(^) largely anti- 
correlates with A2{z). Amazingly, it increases over ~ 
5 Gyr until the peak at ^ ~ 0.8, when l^b starts a sharp 
decay. ft\){z) experiences a break at z ~ 0.35 and be- 
comes shallower, while the bar weakens abruptly (minor 
merger). This evolution is analyzed in §3. 

To view the evolution of stellar populations in the bar 
and the surrounding disk, we have separated the stars 
into age groups: 5-6 Gyrs, 3-4 Gyrs, 2-3 Gyrs, 1-2 Gyrs, 
0-1 Gyrs and 0-0.03 Gyrs old. These are shown via 
frames at z = 0.3 and z = for face-on (Fig. 3) and 
edge-on (Fig. 4) disks. 

Division into populations reveals certain aspects of the 
internal dynamics in the disk. The older stellar popula- 
tion of 5-6 Gyr dominates the disk and bar masses, and 
the ongoing SF is peaked at the center at z = 0. Figs. 3 
and 4 display colors based on stellar ages. The bar has 
a prominent rectangular shape only when observed face- 
on in 1-3 Gyr old stars at z = 0.3. This means that 
the bar extends to the ultra-harmonic resonance and its 
outer part is populated by stars which are trapped by 
the characteristic 4:1 orbits. 

A nuclear bar emerges after z ^ 0.1 (Fig. 3) for 0-3 Gyr 
aged stars at z = 0. One can understand the appearance 
of this morphology because the prime bar slows down 
substantially (Fig. 2) and this leads to an inner Lindblad 
resonance (ILR). A new family of periodic orbits forms 
inside this resonance and these orbits are known to be 
elongated normal to the orbital family which supports 
the main bar figure. This new family of orbits is popu- 
lated by stars with the lowest dispersion velocities, i.e., 
the youngest stars inside the ILR because the trapping 
ability of the nuclear bar is low. 

3. DISCUSSION: DYNAMICS AND POPULATIONS 

We have studied the evolution of galactic bars in fully 
self- consistent cosmological simulations of galaxy forma- 
tion. In a representative model, we find that the first 
generation of bars forms early and in response to the 
asymmetric background DM distribution. This bar de- 
cays quickly and subsequent bars form and are destroyed 
by interactions with subhalos. Finally, a long-lived bar 
is triggered that survives for ~ 10 Gyr. Most interest- 
ingly, its pattern speed increases substantially over the 
first ~ 5 Gyr, despite the angular momentum loss to the 
outer disk and the cuspy DM halo. The bar evolution 
appears to be closely linked to the disk evolution — it 
is permeated by interactions with subhalos as we discuss 
below. 

The sequence of bars in the model are repeatedly trig- 
gered by interactions with DM substructure and do not 
form as a result of a 'classical' (i.e., intrinsic) bar instabil- 
ity. The number of subhalos within the inner prime halo 
is highly variable — the subhalos cluster in the filament 
and are frequently accreted in groups before they merge 
(Romano-Diaz et al. 2008b). This amplifies the dam- 
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Fig. 1. — Penetrating minor merger triggering gas-rich stellar bar in the growing disk. The intruder is on the prograde trajectory inclined 
to the disk plane by about 50° with an impact parameter of ~ 9 kpc. Left: gas; Right: stars — the color palette traces the stellar ages 
(0-0.2 Gyr blue; 0.2-2 Gyr aquamarine). The intruder is also gas rich and is shown near the closest first approach. 
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Fig. 2. — Evolution of the bar amplitude A2 (upper) and pattern 
speed (lower) with the cosmological expansion factor a and 
redshift z. The gap corresponds to a particularly intensive merger 
activity when both curves are ill-defined. is normalized by 

the angular velocity of a test particle at the radius of a maximal 
rotation velocity in the halo at 2; = 0. The first generation bar 
(during the first Gyr) is induced by the main DM filament, while 
the next generations are tidally-induced and destroyed. Finally, a 
long-lived bar is tidally triggered by the interaction with a subhalo 
in a a highly-inclined prograde encounter (Fig. 1). 



age inflicted by the substructure on the disk-bar system 
because of higher mass and energy influx into the disk re- 
gion. We flnd that the prograde encounters trigger bars, 
direct hits weaken existing bars, while retrograde encoun- 
ters have httle effect — in hue with Gerin et al. (1990) 
and Berentzen et al. (2004). In spite of this, addition of 
the cold gas to the disk, up to ~ 15 kpc, in tandem 
with growing stellar surface density by ~ 10, allow the 
bar to speed up over ~ 5 Gyr (Fig. 2). This happens be- 
cause the buildup of a disk mass and the increased central 
mass concentration in the galaxy amplifles the precession 
rate of bar orbits. On the other hand, the weakening of 
the bar during this time results from the angular mo- 



mentum added by the gas to the disk and, therefore, to 
the bar. This happens because orbits trapped within the 
bar become less eccentric, leading to a weaker bar. 

The gas fraction in the disk exceeds that considered so 
far in the literature by as much as a factor of 10 (Bour- 
naud et al. 2005; Berentzen et al. 2007). It is an influx 
of subhalos into the center which strips the disk of its 
gas and heats up the stellar 'fluid' injecting it above the 
plane. Being cut off from the source of angular momen- 
tum, the bar brakes rapidly. 

An important caveat of this slowdown is the bar weak- 
ening after z ^ 0.1 (Fig. 2). Here we do not observe the 
anticorrelation between A2 and The appearance of 
the nuclear bar sheds light on this behavior. The orbit 
trapping by the nuclear bar which grows due to a slow- 
down deprives the main bar of orbital support — the 
weakening of the main bar and the strengthening of the 
nuclear bar proceed in tandem. 

Stellar populations provide a further insight into the 
dynamics of the bar-disk systems. When viewed in col- 
ors of a younger stellar population, the disk is thinner 
and shows less contribution to the spheroidal component. 
Comparison of the left frames in Fig. 4 also shows that 
the spheroidal component is populated by stars older 
than the average disk population, as it disappears in the 
lower frame. 

An interesting feature is present in z = frames of 
Fig. 4 with stars younger than 4 Gyr in the center — 
the upper/lower surfaces are flat and parallel to the disk 
midplane. It is associated with the vertical oscillations 
of stars in the region, and closely resembles the boxy 
bulges observed in about 50% of edge-on disk galax- 
ies (e.g., Lutticke et al. 2000). These bulges originate 
in trapping by a family of 3-D orbits belonging to the 
prime bar (e.g.. Combes et al. 1990; Skokos et al. 2002; 
Martinez- Valpuesta et al. 2006). 

The SF continues after the major merger epoch, z ^ 
1.5 (Fig. 5). The asymptotic value for the SF rate, SFR~ 
1 M0 yr~^, in the disk was not flne-tuned but results 
from self-regulation. Early SF is supported by the cold 
gas accretion onto the disk. In the later stage, the SF 
is limited to the central kpc only — this explains the 
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Fig. 3. — Stellar populations of the face-on disk divided into age groups at z = 0.3 (upper) and z = Q (lower). From left to right the 
ages are: 5-6 Gyr, 3-4 Gyr, 2-3 Gyr, 1-2 Gyr and 0-1 Gyr. The right frames also include the subgroup of 0-0.03 Gyr old stars. Note 
the appearance of a nuclear bar oriented at 90° with respect to the main bar at z = and visible for stars younger than ~ 3 Gyr (lower 
frames). Frames: 17 kpc on the side. 




Fig. 4. — Same as Fig. 3 but for an edge-on disk. 
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Fig. 5. — SF rates shown within the inner 10 kpc of the disk. 



stellar age gradient in the disk which is positive (Figs. 3, 
4). The bar resembles those in the early type galaxies 
with a constant surface brightness along the bar major 
axis, unlike the exponential bars which show negative age 
gradients (Perez et al. 2007). 

In summary, we have followed the bar evolution in the 
most realistic environment attempted so far. We have 
shown that bars can maintain their fast rotation for ^ 
5 Gyr despite the loss of angular momentum to the disk 
and cuspy DM halo. Finally, bar dynamics appears to be 
closely linked to the disk stellar populations. 

We are grateful to our colleagues, and especially to 
John Dubinski, Shardha Jogee, Seppo Laine, Irina Mari- 
nova, Reynier Peletier and Jorge Villa- Vargas for discus- 
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